AD-A122  291  CRACK  TIP  PLASTICITV  ASSOCIATED  WITH  CORROSION  ASSISTED 
FATIGUE(U)  SOUTHWEST  RESEARCH  INST  SAN  ANTONIO  TX 
D  L  DAVIDSON  ET  AL.  21  DEC  82  N8B814-7S-C-1828 

F/G  11/6 


1/1 


UNCLASSIFIED 


NL 


SECURITY  CLASSIFICATION  OF  this  RACK  (Whan  Data  Entarad) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


4.  TITLE  (and  Submit) 

CRACK  TIP  PLASTICITY  ASSOCIATED 
WITH  CORROSION  ASSISTED  FATIGUE 


7.  AU  THORf»> 

D.  L.  Davidson 
J.  Lankford 


9.  performing  organization  name  ano  aooress 

Southwest  Research  Institute 
6220  Culebra  Road,  P.0.  Drawer  28510 


kfirWiTTi  .i.i.u.i  swim:' 


It.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 


S.  TYPE  OF  REPORT  ft  PERIOD  COVERED 

FINAL  REPORT 

June  1974-December  1982 


Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  VA  22217 


S.  CONTRACT  OR  GRANT  NUMBER/*) 


N00014-75-C-1038 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMEERS 


NR  036-109/2-25/76(471) 


12.  REPORT  OATB 


mvsBmiErmiimvm 


.  MONITORING  AGENCY  name  ft  AOORESS/!/  dtttarant  trom  Controlling  Ottlca) 

IS.  SECURITY  CLASS,  (ot  thla  report ) 

| 

( 

ISa.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

IS.  DISTRIBUTION  STATEMENT  (ol  tbit  Report; 

Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the 
United  States  Government.  Distribution  is  unlimited. 


\C 


19.  KEY  WOROS  (Conttnua  on  tararaa  t Ida  it  nacaaaary  and  :damity  by  block  numbar) 


Corrosion  fatigue  Crack  tip  strains  Electron  channeling 

Crack  tip  plasticity  Aluminum  alloys  Subcell 

Fatigue-environment  interaction  Crack  opening  displacement  Crack  tip  stress 
Fractography  Low-carbon  steel 

Fatigue  crack  propagation _ Stereoimaging _ 

20.  ABSTRACT  r Continue  on  ravaraa  net  it  nacaaaary  and  idamlty  by  block  numbar) 

- Measurement  of  the  influence  of  water  vapor  in  the  environment  on  fatigue 

crack  tip  strains  and  stresses  has  been  the  objective  of  this  program.  In  low- 
carbon  steel,  stresses  in  the  plastic  zone  were  deduced  by  imaging  subcells 
formed  during  fatigue  crack  growth.  Water  vapor  was  found  to  lower  the  stress 
levels  and  decrease  the  size  of  the  zone  of  stresses  exceeding  the  yield  value. 
Strains  within  the  plastic  zone  and  at  the  crack  tip  may  be  determined  by  the 
stereoimaging  technique.  For  low-carbon  steel,  water  vapor  lowers  the  strains 
carried  by  the  crack  tip,  and  stresses  computed  from  these  strains  agree  well\^ 


»CRM 
JAN  ?7 


SSlTtCN  OF  l  NOV  49  1  3S30LX~t 


_ UNCLASSIFIED _ 

P'  C  Z  *  "  1 5  *  #n  Z&tm  £“r  » 


UNCLASSIFIED 


SCCUftlTY  CLASSIFICATION  OF  THIS  FAOC(1*h«n  Data  Bn  farad) 


>with  those  determined  from  measurement  of  subcell  size.  A  manuscript  Is 
Included  In  this  report  which  describes  the  computational  technique  of  con¬ 
verting  strains  to  stresses. 
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Strains  In  the  crack  tip  region  have  been  extensively  determined  for 
aluminum  alloys:  7075-T651  ingot  alloy  and  an  experimental  powder  metallurgy 
alloy  MA-87.  Water  vapor  in  the  environment  is  shown  to  have  a  similar  effect 
on  strains  as  in  low-carbon  steel.  Dynamic  observations  of  crack  growth  in 
vacuum  have  revealed  that  cracks  do  not  grow  on  each  cycle.  Fractography 
supports  both  of  these  observations,  and  reveals  many  other  details  of  crack 
growth  as  well. 

Detailed  crack  tip  measurements  have  revealed  the  fact  that  crack  opening 
is  mixed  mode  even  though  loading  is  purely  Mode  I.  The  degree  of  Mode  II  in 
the  crack  opening  is  strongly  influenced  by  the  environment,  level  of  cyclic 
stress  intensity  factor,  and  the  metallurgical  state  of  the  materials.  This 
finding  is  potentially  a  very  important  outcome  of  this  research. 

A  significant  effort  was  expended  to  utilize  the  factors  measured  in 
models  for  crack  growth,  and  some  progress  was  made. 
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SCOPE  AND  PURPOSE  OF  INVESTIGATION 


Fatigue  crack  growth  has  been  shown  to  be  dependent  on  environment, 

especially  at  values  of  cyclic  stress  intensity  factor  (AK)  approaching 

r  r  y* 

the  threshold  value.  '  This  behavior  has  been  termed  "environmentally- 
assisted  fatigue  crack  growth" to  separate  this  behavior  from  that 
related  to  cyclic  stress  corrosion.  This  program  was  planned  to  utilize 
new  tools  for  measuring  the  effect  of  environment  on  the  material  response 
(plasticity)  associated  with  fatigue  crack  growth. 

Measurements  were  made  on  the  effect  of  environment  and  AK  on  crack 
tip  plasticity  parameters  for  low  carbon  steel ,(A»B»C»D»1 »2,3)  7075.76 
ingot  aluminum  alloy, (F>G»10»13,17)  ^-87  a  powder  metallurgy  derived 
aluminum  alloy, »13,16,17)  an<j  in  the  recrystalllzed- 

annealed  condition. 

Three  environments  were  used:  very  dry  nitrogen  (less  than  10  ppm 
water  vapor),  vacuum  of  10“8  torr  (1.3  mPa)  and  humid  air  (not  less  than 
12,000  ppm  water  vapor).  Differences  between  vacuum  and  very  dry  nitrogen 
results  could  not  be  detected. 

Two  techniques  were  used  to  study  crack  tip  plasticity.  At  the 
start  of  the  investigation,  low-carbon  steel  was  found  to  form  subcells 
near  fatigue  cracks,^*8,1  *2^  and  these  were  readily  observed  by  electron 
channeling  contrast  In  the  SEM.^®*28^  This  technique  was  used  to 
characterize  the  subcell  distribution  as  a  function  of  environment  and 
AK,  and  the  Information  was  used  to  derive  both  the  stress  distribution 
near  the  crack  tip,  and  to  compute  the  energy  dissipated  during  crack 
growth. (B»c»5»6»7»8)  j0Wards  the  end  of  this  part  of  the  Investigation, 
two  new  techniques  were  developed  for  studying  fatigue  cracks:  a  cyclic 
loading  stage  was  developed  for  the  SEM^^1^  so  that  high-resolution, 
real-time  observations  of  crack  growth  could  be  made,  and  the  stereoimaging 
technique ^ IV^  was  developed  for  accurate,  high-resolution  determination  of 
crack  tip  strains  and  opening  displacements.  These  new  techniques  were 
then  applied  to  all  the  materials  investigated. 

* 

A-G  refer  to  technical  reports;  journal  references  reporting  work  from 
this  contract  are  numbered  1-21;  Roman  numerals  are  used  to  designate 
journal  references  not  specifically  related  to  this  contract. 
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SUMMARY  OF  RESULTS 

The  general  finding  of  this  investigation  common  to  all  the 
materials  studied  is  that  water  vapor  In  the  environment  decreases  the 
strain  range  which  the  crack  tip  can  support, (7,10,11,14)  though  for 
the  titanium  alloy,  this  finding  is  subject  to  interpretation.  A  typical 
result  for  7075-T651  is  shown  In  Fig.  1.  Crack  tip  strain  and  crack 
opening  are  found  to  be  strongly  dependent  on  AK  for  both  environments 
and  all  materials.  For  near- threshold  AK,  crack  tip  strain  approaches 
the  elastic  limit  in  the  absence  of  water  vapor,  and  there  is  some  evi¬ 
dence,  although  not  conclusive,  that  crack  growth  may  occur  at  strains 
below  the  elastic  limit  in  the  water  vapor  environment.  Again,  Ti-6A1-4V 
does  not  completely  fit  this  pattern.  Another  Important  finding  of  this 
research  is  that  cracks  do  not  grow  on  each  cycle  in  the  near-threshold 
region/10*11  ^  although  in  humid  air,  growth  occurs  more  frequently, 
i.e.,  after  fewer  number  of  cycles,  than  In  the  absence  of  water  vapor. 

Crack  opening  has  been  found  to  be  mixed  mode  in  each  alloy  system 
investigated  and  under  all  conditions  of  stress  intensity  and  environ¬ 
ment/17^  But  the  magnitude  of  each  mode  is  strongly  influenced  by  each 
of  these  variables.  Models  of  crack  tip  plasticity  and  opening  do  not 
presently  consider  Mode  II  (shear)  as  being  a  significant  factor,  yet  in 
MA-87  this  is  the  dominant  opening  mode/17^  The  influence  of  metal¬ 
lurgical  structure  on  this  aspect  of  crack  tip  mechanics  Is  potentially 
an  important  factor  in  controlling  fatigue  crack  growth,  and  needs  to  be 
further  explored. 

The  energy  dissipated  during  fatigue  crack  growth  was  computed  for 
low  carbon  stee/8*0,5’8’7’8)  and  found  to  depend  on  environment  and 
stress  Intensity.  These  observations  and  computations  were  incorporated 
Into  a  model  for  near- threshold  crack  growth/8^  Water  vapor  in  the 
environment  results  in  a  lower  rate  of  energy  dissipation  per  unit  of 
crack  formation. 

Models  for  crack  growth  in  the  aluminum  and  titanium  alloys  have 
been  developed^8*10*11 ’  ^  using  information  determined  by  the  stereo- 
imaging  technique. 


Met  Environment 


Distribution  of  maximum  shear  strain  range  for  7075-T651  in  both 
wet  and  dry  environments  at  AK  *  6  MN/m3/2.  Crack  is  shown  on 
the  plane  of  zero  strain.  X  and  y  are  in  micrometers.  Scale  on 
strain  axis  is  the  same  for  both  plots.  Water  vapor  in  the 
environment  causes  a  reduction  in  the  strain  range. 
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The  results  discussed  above  are  only  the  most  general  conclusions 
reached;  there  were  numerous  interesting  experimental  findings  related  to 
both  fatigue  crack  growth,  the  effect  of  environment,  and  the  mathematical 
modeling  of  fatigue  crack  growth.  These  results  have  been  published  in 
the  form  of  annual  technical  reports  and  papers.  The  technical  reports 
have  also  been  used  to  disseminate  more  detailed  experimental  results  than 
it  has  been  possible  to  put  into  journal  articles.  The  published  papers 
are  more  complete  and  interpretive  than  the  reports.  Between  these  two 
forms  of  communication,  nearly  all  the  results  derived  during  this  research 
have  been  reported. 


PROJECT  SUMMARY 

Two  new  techniques  for  the  study  of  environmental  alteration  of 
material  plasticity  at  crack  tips  have  been  developed  during  the  life  of 
this  contract.  These  techniques  allow  the  determination  of  material 
stresses  and  strains  within  a  very  small  region  close  to  the  crack  tip 
(1-10  ym).  These  measurement  techniques  were  used  to  determine  the  effect 
of  environment  on  fatigue  crack  growth  in  low  carbon  steel,  the  aluminum 
alloys  7075-T651  and  MA-87  (a  P/M  product)  and  Ti-6A1-4V  (RA).  The  in¬ 
formation  gained  by  these  measurements  was  used  in  models  whose  objectives 
were  to  both  correlate  the  various  parameters  measured  and  to  gain  more 
fundamental  insight  into  the  factors  controlling  environmentally-assisted 
fatigue  crack  growth. 


CONCLUSIONS 

The  strain  which  a  crack  tip  can  support  is  lowered  when  the  crack 
is  grown  in  a  water  vapor  environment,  as  compared  to  a  very  dry  or  vacuum 
environment.  Hydrogen  formed  by  dissociation  of  water  vapor  at  the  newly 
formed  surface ^ is  believed  to  cause  this  "embrittlement."  All  indica¬ 
tions  are  that  the  hydrogen  is  acting  within  10  ym  of  the  crack  tip,  or 


VV.vVV.V-> 


less.  In  aluminum  alloys,  the  details  of  this  effect  are  influenced  by 
heat  treatment  and  microstructure. 
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TABLE  I 

CRACK  TIP  OPENING  PARAMv  ERS  1  urn  BEHIND  CRACK  TIP, 
AK  3  7.5  MNm-3/2 


CTOD* 


7075-T6 

7075-T6 

MA-87 


MA-87 


Environment 


Vacuum 


Vacuum 


0.238 
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The  fractography  of  7075-T651  is  conventional  over  this  stress 

5 

intensity  regime,  i.e.,  striations  are  formed  in  both  air  and  vacuum. 

The  behavior  of  MA-87  is  markedly  different.  In  vacuum,  for  example, 
fatigue  crack  growth  produces  approximately  spherical  particles,  which  are 
dispersed  over  a  fractographic  background  of  diffuse,  ripply  striations; 
these  features  are  observed  uniformly  across  the  specimen  fracture  surface. 
The  particles  usually  are  observed  as  Isolated  individuals,  but  they  fre¬ 
quently  occur  in  clusters  (Figures  2  and  3)  as  well.  Energy  dispersive 
X-ray  analysis  (EDS)  of  several  of  the  spheres  showed  that  they  have  the 
same  chemical  makeup  as  the  matrix  material.  Thus,  the  spheres  do  not 
correspond  to  the  small  (0.02  pm  to  0.3  \m)  ^Alg  dispersoid  particles6 


characteristic  of  MA-87.  However,  It  Is  possible  to  derive  some  idea  as 
to  their  origin  by  considering  clusters  of  the  spheres  in  various  stages 
of  development. 

For  example.  In  Figure  2  can  be  seen  a  large,  smooth  platelet  (P), 
from  which  several  connected  spheres  (A)  have  separated.  Since  it  appears 
that  the  platelet  and  the  spheres  once  fit  together,  the  spheres  evidently 
formed  while  still  attached  to  the  platelet.  Other  spheres,  such  as  B, 
show  what  appear  to  be  laps,  as  though  they  had  been  flattened  out  and  then 
rolled  up.  Once  a  platelet  Is  formed,  parts  of  it  evidently  begin  to  roll 
up;  these  eventually  detach  from  the  parent  platelet  by  fracture,  to  form 
separate  spheres  such  as  C  in  Figure  3.  Note  the  apparent  rolled-up  pre¬ 
sphere  on  platelet  P  in  this  figure. 

It  would  appear  that  the  critical  step  In  sphere  formation  is  the 
creation  of  the  platelets.  Such  features  have  in  fact  been  observed 
before,  but  only  in  cases  in  which  Mode  II  interfacial  sliding  was  ex¬ 
ternally  imposed  by  the  specific  loading  conditions.  These  have  included 
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fretting  fatigue,  rolling  contact,  and  pure  Mode  II  fatigue,  i.e., 

cyclic  compression  of  slant  cracks.  Under  such  circumstances,  localized 

regions  of  opposing  surfaces  experience  sliding  contact,  each  producing 

an  extraordinarily  smooth,  thin  layer  of  material  known  as  a  Beilby 

layer.™  The  Beilby  layer  has  the  same  chemical  composition  as  the  parent 

alloy,  but  crystalline  order  Is  almost  destroyed  by  the  rubbing  process, 

and  the  end  product  Is  extremely  ductile.  It  is  this  anomalous  ductility 

which  permits  the  layers  to  roll  up  into  laminated  spheres.™ 


It  is  Interesting  that  spheroidal  debris  Is  not  found  on  MA-87 
fatigue  fractures  formed  In  air.  Instead,  the  surface  Is  heavily  oxidized 
(Figure  4),  so  much  so  that  strlatlons,  although  present,  are  hard  to 
discern.  Superimposed  on  the  fracture  surface  are  very  small,  round  par¬ 
ticles,  seemingly  embedded  In  the  fracture  surface.  These  are  much  smaller 
than  the  spherical  wear  particles  formed  in  vacuum,  and  are  of  the  size, 
shape,  and  spacing  of  the  Co2Alg  dispersoids.  Electron  beam  diameter  and 
penetration  into  the  underlying  substrate  preclude  EDS  determination  of 
their  composition. 

Irregular-shaped  bits  of  apparent  oxide  debris  (arrow.  Figure  4) 
also  are  evident.  Their  presence  Is  best  appreciated  by  examining  the 
flanks  of  cracks  grown  In  air;  these  are  seen  to  exude  debris  (Figure  5). 
The  electron-charging  nature  of  the  debris  strongly  suggests  that  it 
derives  from  fractured  oxide  films. 

By  considering  certain  recent,  related  research  involving  similar 
alloys,  it  is  possible  to  understand  the  tendency  of  MA-87,  but  not  7075- 
T651,  to  produce  wear  debris,  and  for  the  MA-87  debris  particles  to  be 
smooth  and  spherical  In  vacuum,  while  rough  and  Irregular  In  air.  In  the 
first  case,  Vasudevan  and  Suresh^  have  recently  Investigated  the  oxide- 
forming  propensity  of  7XXX  aluminum  alloys,  of  which  MA-87  is  a  variant. 
Among  their  findings  is  the  fact  that  during  fatigue  crack  growth  in  air, 
the  level  of  oxide  formation  on  fresh  crack  surface  is  on  the  order  of 
0.01  pm,  negligible  compared  to  crack  tip  opening  displacements  (Table  I), 
for  peak-aged  alloys  like  7075-T651.  However,  for  overaged  alloys  like 


MA-87*,  oxide  layers  of  up  to  0.1  pm  In  thickness  are  produced  during 

fatigue  crack  growth  in  moist  air  at  near- threshold  stress  intensities. 
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Rabinowicz  has  noted  that  oxidation  during  sliding  contact  will 
generally  lead  to  the  formation  of  debris  which  prevents  two-body  burnish¬ 
ing.  The  debris  occurs  by  the  fracture  and  spalling-off  of  the  oxide 
layer.  In  the  present  case,  this  process  apparently  accounts  for  the 
debris  which  is  extruded  during  the  growth  of  cracks  in  MA-87  in  air;  it 
will  be  recalled  (Table  I)  that  the  Mode  II  crack  tip  displacement  com¬ 
ponent  is  quite  large  relative  to  the  total  CTOD.  For  7075-T651,  the 
oxide  layer  formed  in  air  is  negligible,  and  so  is  the  Mode  II  opening 
displacement;  hence,  no  oxide  debris  is  observed. 

On  the  other  hand,  there  is  no  oxide  layer  to  prevent  the  local 
burnishing  of  MA-87  crack  interfaces  in  vacuum,  and  together  with  the 
large  Mode  II  crack  tip  displacement,  the  formation  of  Beil  by  layers  and 
their  transformation  to  spheroidal  debris  is  favored.  Although  there  is 
also  no  oxide  to  prevent  burnishing  of  7075-T651  in  vacuum,  the  Mode  II 
displacement  in  this  case  is  insufficient  to  produce  the  requisite  damage 
layers.  Comparing  the  CTODjj  values  for  7075-T651  and  MA-87  in  vacuum 
(Table  I),  it  appears  that  the  threshold  Mode  II  displacement  for  the 
production  of  spherical  wear  particles  in  7XXX  aluminum  alloys  must  lie 
between  0.1  pm  and  0.72  ym.  It  also  is  interesting  to  note  that  the 
largest  MA-87  spheroidal  particles  are  on  the  same  order  of  size  as  CTODjj. 

*Heat  treatment  was  -T6  plus  14  hours  at  162C. 


If  these  rationalizations  are  correct,  then  It  must  also  be  true 
that  the  crack  faces  do  not  open  In  Mode  I  before  at  least  some  sliding 
has  occurred.  To  verify  this,  the  opening  of  MA-87  crack  tips  was  photo¬ 
graphed  at  various  Increments  In  the  load  cycle,  and  the  resulting  micro¬ 
graphs  were  analyzed  using  stereoimaging.  It  was  found  that  as  the  cyclic 
stress  Increases  from  Its  minimum  value,  the  crack  tip  starts  sliding  at 
a  very  low  stress  Intensity;  Mode  1  opening  Is  not  achieved  until 
AK  =  0.5  K|nax.  Thus,  there  Is  ample  opportunity  for  Interfaclal  contact 
during  Mode  II  sliding. 

At  present.  It  Is  not  at  all  clear  why.  fatigue  cracks  In  MA-87 
choose  to  open  In  mixed  Modes  I  and  II,  while  7075-T651  cracks  open  pre¬ 
dominantly  In  Mode  I;-  neither  Is  It  apparent  why  the  CTOD  for  MA-87  Is  so 
much  greater  than  that  for  7075-T651.  Aside  from  grain  size,  ~50  vim  for 
7075-T651  versus  -5  ym  for  MA-87,  the  principal  differences  between  the 
two  alloys  lie  In  their  microstructures.  MA-87  Is  overaged,  and  has  a 
low  dlspersold  population,  while  7075-T651  Is  peak-aged,  with  a  high  dls- 
persold  content. 

Regardless  of  the  origin  of  the  Mode  II  opening  In  MA-87,  Its 
presence  Is  remarkable,  considering  the  relatively  low  driving  force 
responsible  for  It.  Analysis  of  the  present  crack  profiles  shows  that  on 
average,  for  both  alloys  In  air  and  vacuum,  a  given  crack  segment  Is 
oriented  approximately  25°  from  normal  to  the  stress  axis.  Under  such 
circumstances,^  the  Mode  I  stress  Intensity  normal  to  the  local  crack 
plane  Is  .82  of  the  Mode  I  stress  Intensity  (Kj)  for  the  average  crack 
plane.  Correspondingly,  the  local  Mode  II  stress  intensity  Is  only  .38  Kj. 


The  presence  of  the  wear  debris  itself  has  implications  for  fatigue 
crack  growth,  since  Vasudevan  and  Suresh^  have  shown  that  the  presence  of 
oxide  layers  alone  can  significantly  alter  the  effective  Mode  I  stress 
intensity  at  low  rates  of  growth.  The  existence  of  the  oxide  wear  debris 
(Figure  5)  means  that  the  wedging  effect  of  the  oxide  layer  is  further  en¬ 
hanced  in  materials  which  have  a  large  Mode  II  crack  opening  component, 
since  the  fracturing  of  the  oxide  layers  will  promote  the  formation  of 
new  oxide,  effectively  increasing  the  thickness  of  the  layer.  Similarly, 
the  presence  of  the  spherical  debris  in  vacuum  corresponds  to  a  wedging 
influence  which  will  likewise  tend  to  reduce  the  Mode  I  effective  stress 
intensity. 
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These  conclusions  are  supported  by  the  work  of  Ryder,  et  al,  who 
studied  the  behavior  of  torsional  Stage  I  cracks  in  an  Al-Mg-Zn  alloy. 

They  found  that  Stage  I  (Mode  II)  crack  growth  rates  were  enhanced  by 
inert  liquids  which  normally  retard  Stage  II  (Mode  I)  crack  growth.  The 
effect  of  the  liquids  was  to  pump  out  of  the  crack  fretting  debris  pro¬ 
duced  by  Mode  II  sliding,  thereby  Increasing  the  effective  stress  intensity. 
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Figure  2.  Spheroidal  wear  debris  on  fatigue  fracture  surface  of 
MA-87,  vacuum  environment,  AK  =  7  MNm-3/2.  Platelet 
(P),  separated  spheroids  (A),  lapped  spheroid  (B). 
Large  arrow  shows  direction  of  crack  growth. 


Figure  3.  Spheroidal  wear  debris  on  fatigue  fracture  surface  of 
MA-87,  vacuum  environment,  AK  =  9  MNm-3/2.  Platelet 
(P)  with  rolled  up  pre-sphere,  spheroid  (C)  detached 
from  matrix.  Large  arrow  shows  direction  of  crack 


Fatigue  fracture  surface  of  MA-87  cycled  at 
AK  =  8  MNm"3/2  jn  air,  showing  heavily  oxidized 
striations.  Large  arrow  shows  direction  of 
crack  growth;  small  arrow  denotes  typical  oxide 
particle. 


Oxide  debris  extruded  from  crack  flank 
MA-87  cycled  in  air,  AK  -  10  MNm'3/2. 
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CRACK  TIP  STRESSES  AS  COMPUTED  FROM  STRAINS 
,  DETERMINED  BY  STEREOIMAGING 


P.  L.  Davidson*,  D.  R.  Will lams +  and  J.  E.  Buckingham* 


Measurements  of  the  tte'ee  in-plane  elements  of  the  strain  tensor  near 
a  crack  ttp  have  teen  combined  with  the  material  constitutive  behavior  to 
compute  the  three  In-plane  elements  of  the  stress  tensor  and  the  mean 
stress,  the  method  is  Illustrated  with  fatigue  crack  plasticity  results 
for  7075-T6  aluminum  alloy. 

INTRODUCTION 

too  of  the  reasons  a  crock  presents  a  threat  to  an  engineering 

is  because  the  nominal  (design)  stress  to  which  the  body  Is  sub¬ 
jected  is  raised  at  the  crack  tip;  this  local  stress  field  is  lit  turn  the 
driving  force  for  crack  extension.  Thus,  a  large  amount  of  effort  has  been 
expended  over  a  long  period  of  time  to  determine  the  magnitude  and  extent 
of  this  stress  concentration.  The  region  of  elevated  stress  at  a  crack  tip 
is  very  small  and  varies  rapidly,  making  experimental  determination  of  Its 
magnitude  a  difficult  task.  It  IS  not  surprising,  therafore,  that  most  of 
the  effort  has  been  theoretical  and  computational.  Fatigue  cracks  present 
special  problems  because  of  the  necessity  of  knowing  the  constitutive 
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behavior  of  materials  subjected  to  cyclic  loading;  therefore,  most  of  the 
analytical  effort  has  gone  into  monotonically  loaded  cracks.  Thus,  it  is 
understandable  that  fatigue  crack  tip  stress  remains  a  poorly  quantified 
parameter.  This  paper  presents  a  method  for  computing  stresses  in  the 
near  field  of  a  fatigue  crack  from  experimentally  determined  strains. 

EXPERIMENTAL  DETAILS 

Single  edge  notched  specimens  of  7075-T6  aluminum  alloy  having  a 

gauge  section  of  20  mm  x  3  mm  thick  were  pin  loaded  to  obtain  crack  tip 
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stress  intensities  of  6  <  AK  <  11  MN/m  '  at  a  frequency  of  1  Hz  in  a 
vacuum  of  1.3  mPa.  The  direction  of  crack  growth  was  perpendicular  to 
the  rolling  direction,  and  crack  lengths  5  <  l  <  10  mm  were  investigated. 

DETERMINATION  OF  STRAINS 

Visualization  and  measurement  of  crack  tip  displacements  with  high 
spatial  resolution  are  possible  using  the  stereoimaging  technique. ^ 

This  technique  employs  the  ability  of  the  eye-brain  system  to  detect  dif¬ 
ferences  between  two  views  of  an  object  by  constructing  a  three-dimensional 
image  of  the  object.  For  crack  tips,  it  is  only  necessary  to  obtain 
photographs  in  two  loading  states  in  order  to  image  the  displacements 
which  occurred  during  the  loading  increment.  The  displacements  imaged 
are  in  the  plane  of  the  photograph  and  along  the  eye  axis.  These  displace¬ 
ments  can  be  quantified  by  using  the  techniques  of  photogrammetry.  The 

measured  displacements  then  may  be  used  to  calculate  three  elements  of  the 

(2) 

symmetric  in-plane  strain  tensor.'  '  Very  high  spatial  resolution  around 


the  crack  tip  can  be  obtained  by  using  a  special  cyclic  loading  stage 
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incorporated  into  the  specimen  chamber  of  a  scanning  electron  microscope.  ' 
In  this  manner,  the  high  resolution  and  depth  of  field  characteristics  of 
the  SEM  are  utilized.  Strains  of  up  to  nearly  100%  have  been  measured, 
within  1  uni  of  the  crack  tip  using  the  SEM  fatigue  stage  and  stereoimaging 
technique.  For  this  purpose,  photographs  at  magnifications  up  to  about 
10,000X  have  been  made,  although  such  high  magnification  is  not  normally 
usable  because  of  the  limited  area  covered.  A  pair  of  typical  crack  tip 
photographs  is  shown  in  Figure  1.  This  pair  of  photographs  was  used  to 
determine  the  displacements  shown  in  Figure  2.  The  asymmetry  of  dis¬ 
placements  shown  in  this  figure  is  typical  of  fatigue  cracks  grown  at 
near- threshold  stress  intensities,  and  it  is  believed  to  be  due  to  the 
propensity  of  metals  to  deform  in  shear.  At  higher  stress  intensities, 
greater  symmetry  of  displacements  is  observed.  As  shown  in  Figure  la,  the 
crack  is  closed  at  minimum  load;  crack  opening  for  this  case  has  been 
found  to  occur  at  approximately  0.4  of  maximum  load.  At  full  load,  crack 
opening  is  a  mix  of  Mode  I  and  Mode  II,  as  may  be  deduced  from  displace¬ 
ments  behind  the  crack  tip  in  Figure  2. 

From  the  x  and  y  displacements  of  Figure  2,  the  coordinate  strain 

increments*  AX  and  AX  ,  and  the  shear  strain  increment  Ay  may  be 
xx  yy  xy 

computed. ^  The  principal  strain  increments  AX-j  and  AX 2  and  the  maximum 
shear  strain  increment  Aymax  may  then  be  computed  from  the  coordinate 

* 

The  strains  computed  are  referred  to  as  strain  increments  because  they 
refer  to  the  strain  change  between  two  loading  states.  For  a  fatigue 
crack  tip,  the  strain  at  minimum  load  is  usually  not  at  zero  strain, 
as  referenced  to  the  uncracked  state. 


strains  by  the  appropriate  transformation  operations  (Mohr's  circle). 

These  in-plane  strains  are  all  that  may  be  determined,  since  no  method  has 

yet  been  devised  to  measure  the  z-axis  displacement  gradients.  Typical 

distributions  of  the  principal  strain  increments  are  shown  in  Figures  3a, 

b,  and  c.  These  were  derived  from  the  displacements  of  Figure  2.  A 
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detailed  account  of  computation  of  strains  has  already  been  given.  ‘ 

COMPUTATION  OF  STRESSES 

The  procedure  for  computing  stresses  from  strains  was  derived  by 
(41 

Lubahn  and  Felgar.  '  Their  procedure  was  developed  to  use  strains  as 
measured  by  electrical  resistance  strain  gage  rosettes  placed  on  the  sur¬ 
face.  The  use  of  their  procedure  with  strains  determined  from  stereo¬ 
imaging  is,  therefore,  a  completely  analogous  situation. 

The  surface  of  any  specimen  is,  by  definition,  in  a  state  of  plane 
stress,  which  means  that  no  constraint  exists  normal  to  the  surface.  Thus, 
a  limitation  of  the  stereoimaging  technique  is  that  only  surface  stresses 
may  be  computed. 

The  strains  determined  by  stereoimaging  are  total  strains,  and,  since 
they  are  the  result  of  changes  between  two  loading  states,  they  are  actually 
total  strain  increments;  thus,  L\  • ,  the  measured  strain  increment,  is  in 
general 

A*.  =  Ae.  +  Aei  (1) 

where  Ae^  and  Aei  are  the  plastic  and  elastic  strain  increments,  respectively. 
With  several  assumptions,  Appendix  A,  it  is  possible  to  separate  the  elastic 
and  plastic  strain  increments  and  to  compute  the  principal  stress  increments 


Aa-i  *  Aa?  and  at  ,  and  the  coordinate  stress  increments  a ov¥,  a c  and 

I  L  mo  a  xx  yy 

At  .  A  constitutive  equation  is  required  to  make  the  strain-to-stress 
xy 

transformation,  and  it  has  been  necessary  to  assume  that  the  principal 
stresses  have  the  same  direction  as  the  principal  strains  over  one-half  of 
the  loading  cycle. 

The  computational  procedure  was  devised  by  Lubahn  and  Felgar  for 
nonreverse- loaded  bodies,  but  its  use  for  cyclic  loading  is  believed  to  be 
appropriate  for  several  reasons:  (1)  the  strain-to-stress  transformation 
is  made  only  for  each  half  cycle;  therefore,  the  increment  of  straining  is 
for  only  the  tensile  or  compressive  loading  segment  of  the  cycle;  and  (2) 
the  assumptions  of  plasticity  theory  are  equally  valid  in  tension  and 
compression. 

The  strain  increment  data  of  Figure  3  have  been  converted  to 
incremental  stress  using  the  procedure  given  in  Appendix  A,  incorporating 
the  mechanical  properties  values  in  Table  I.  Maximum  principal  stress  in¬ 
crement  (ao-|),  minimum  principal  stress  increment  (a and  the  maximum 
shear  stress  increment  (Axmax)  distributions  near  the  fatigue  crack  tip  are 
shown  in  Figures  4a,  b,  and  c. 

The  mean  stress  is  also  easily  computed  from  Aa-|  and  Aag.  and  the 
distribution  of  this  term  is  shown  in  Figure  5.  Mean  stresses  are  often 
considered  to  be  useful  in  understanding  the  effects  of  environment  and 
stress  state  on  crack  tip  plasticity. 

The  conversion  of  the  principal  stresses  act^  back  into  coordinate 
stresses  in  the  fixed  spatial  x-y  coordinate  system  is  easily  accomplished 
by  using  a  Mohr's  circle  construction.  Coordinate  stresses  are  required  in 


cases  where  a  sequence  of  loading  and  unloading  is  followed  step  by  step, 
in  that  the  summation  of  stresses  may  be  done  only  in  the  fixed  coordinate 
system. 


TABLE  I 

Mechanical  Properties 
7075-T6 

Tensile  Proportional  Limit  Stress  448  MPa 
Tensile  0.2%  Offset  Yield  Stress  508  MPa 

Stress-Strain  Curve  in  Tension  Ao(MPa)  *  7.1xl04  ax  0  <  AX  <  .0166 

a  -062 

A0(MPa)  *  1620(^p)  .0166  -  Ae  <  1 

DISCUSSION 

The  procedure  for  computing  stresses  from  strains  outlined  herein  is 
complex  only  to  the  degree  expressed  by  Equations  (A5)  and  (A6).  The 
quadratic  form  of  these  equations  does  pose  some  problems  with  convergence, 
but  these  have  been  circumvented  by  anticipating  the  conditions  for  which 
convergence  had  been  achieved.  A  procedure  was  formulated  to  force  con¬ 
vergence,  and  it  was  successful  for  the  vast  majority  of  cases  tested. 

Most  of  the  cases  converged  within  3  to  5  iterations  through  the  procedure 
given  in  Appendix  A. 


■  IT  l  .  ,  ■■  _• 


The  significant  opportunities  for  error  in  the  stress  computation 
procedure  are  related  to  the  assumptions  made.  Of  these,  perhaps  the  most 
questionnable  assumption  is  that  of  the  deformation  theory  of  plasticity. 
Equation  (A2),  which  is  valid  only  for  proportional  loading.  Limited  effort 
to  examine  this  issue  indicates  that  strain  actually  does  increase  in  a 
reasonably  proportional  way  for  many  of  the  points  within  the  plastic  zone 
on  the  loading  portion  of  the  cycle. 

Another  possible  source  of  error  is  the  constitutive  equation,  in 
this  case  the  tensile  stress-strain  curve  obtained  on  the  loading  portion 
of  the  cycle  after  a  schedule  of  increasing  strain-increment,  reversed 
cyclic  loading.  The  assumption  is  that  this  curve,  derived  from  uniaxial 
tests,  in  which  aa  is  controlled,  is  applicable  to  the  biaxial  deformation 
of  material  near  the  crack  tip.  Lambda  and  Sidebottom^’^  have  examined 
cyclic  hardening  under  changing  biaxial  conditions,  and  found  that  the 
result  was  essentially  the  same  as  for  uniaxial  cycling  when  compared  on 
effective  stress  and  strain  basis.  Fatigue  crack  movement  in  the  growth 
rate  region  being  studied  (10“®  m/cy  and  below)  is  certainly  slow  enough 
for  a  small  volume  of  material  to  be  cycled  enough  times  to  achieve  a  stable 
a a.  The  other  assumption  related  to  the  tensile  stress-strain  curve  used 
is  that  the  strain  magnitude  of  the  uniaxial  tests  is  lower  than  the  strain 
magnitudes  experienced  near  the  crack  tip,  so  the  measured  results  must  be 
extrapolated  to  large  strains.  The  cyclic  stress -strain  curve,  that  is  the 
curve  relating  the  saturation  value  of  Ao  to  Ae,  as  derived  from  a  strain- 
controlled,  reverse  cycling  uniaxial  specimen,  was  also  used  for  computation 
of  the  stress  for  comparison  with  the  computation  described  in  the  previous 
paragraph.  The  result  produced  values  near  the  crack  tip  which  were  3.4% 

1 ower . 


The  validity  of  using  these  constitutive  equations  near  a  fatigue 
crack  tip  has  not  been  confirmed,  partly  because  there  is  no  independent 
way  to  measure  stresses  near  crack  tips  with  the  required  resolution. 
However,  stresses  computed  by  the  methods  outlined  here  have  been  matched 
with  those  determined  indirectly  from  using  subgrains  formed  near  fatigue 
crack  tips  in  low-carbon  steel  as  "stress  gages,"  and  the  results  were 
found  to  be  very  nearly  the  same.^  Also,  strains  measured  on  a  poly¬ 
carbonate  tensile  specimen  were  converted  to  stresses  using  Hookes  Law  and 
matched  to  the  measured  stress.  For  a  strain  of  approximately  0.01,  the 
stress  computed  was  within  3%  of  that  measured. 

The  results  presented  here  are  interesting,  from  the  viewpoint  of 
crack  tip  mechanics.  The  maximum  principal  stress  increment  at  the  crack 
tip  Is  1490  MPa.  This  is  3  times  the  0.2%  offset  stress  from  a  uniaxial 
tensile  test  (508  MPa),  and  1.7  times  the  cyclic  yield  stress  (880  MPa), 
defined  as  the  stress  for  Ae  =  0.02%.  The  maximum  value  of  mean  stress 
(1180  MPa)  occurred  10  urn  ahead  of  the  crack  tip,  decreasing  to  1120  MPa 
at  the  crack  tip. 

The  outer  plastic  zone,  defined  as  the  distance  ahead  of  the  crack 
tip  where  Acr-|  exceeds  508  MPa,  occurs  at  a  distance  beyond  that  of  the 
strain  analysis,  but  it  is  extrapolated  to  be  110  pm.  The  inner,  or  cyclic 
plastic  zone,  defined  as  the  distance  ahead  of  the  crack  tip  where  Ao-j  ex¬ 
ceeds  880  MPa,  is  60  pm.  These  compare  with  the  net  section  stress  of  50 
MPa  which  would  be  reached  at  a  distance  of  160  pm  ahead  of  the  crack  tip 
(as  determined  by  extrapolation).  Note  that  the  stress  gradient  within 
the  plastic  zone  is  much  more  gradual  than  the  strain  gradient,  as  would 
be  expected  from  the  low  work  hardening  coefficient  of  this  alloy. 


Values  of  plastic  zone  size  (PZS)  may  be  computed  using  the 
equation 

K  2 

PZS  -  <*(£-)  (2) 

a1 

where  PZS  refers  to  the  cyclic  zone  when  =  cyclic  yield  stress  and  to 
the  outer  PZS  when  *  tensile  yield  stress.  This  equation  is  an  esti- 
mate  made  by  Ricev  '  for  non-work  hardening  materials.  He  estimated  that 
for  plane  stress  aQ  =  0.32  for  the  outer  plastic  zone  size  ahead  of  the 
crack,  and  ctc  »  ^  aQ  for  the  cyclic  plastic  zone.  For  the  plastic  zone 
size  estimates  given  in  the  previous  paragraph,  aQ  and  ac  are  computed 
to  be:  aQ  =  .44,  ac  =  .72,  with  ac  *  1.6  aQ.  Estimates  of  PZS  by 
Eq.  (2),  therefore,  result  in  values  smaller  than  are  determined  here. 

SUMMARY  AND  CONCLUSIONS 

1.  A  procedure  for  computing  the  three  in-plane  elements  of  the 

principal  stress  tensor  from  the  in-plane  strain  tensor,  as 

measured  by  stereoimaging,  has  been  adapted  from  the  work  of 

(41 

Lubahn  and  Felgar. N  ' 

2.  Stresses  in  the  near-crack  tip  region  have  been  computed  for 
a  fatigue  crack  tip  in  7075-T6.  The  maximum  and  minimum 
principal  stresses,  maximum  shear  stress  and  mean  stress  dis¬ 
tributions  near  the  crack  tip  are  illustrated,  together  with 
the  principal  strain  distributions,  which  are  the  measured 
values. 

3.  The  maximum  principal  stress  is  computed  to  be  3  times  the 
yield  stress  and  1.7  times  the  cyclic  yield  stress.  The 
maximum  value  of  mean  stress  is  found  to  be  10  urn  ahead  of 
the  crack  tip. 


APPENDIX  A 


Calculation  of  the  principal  stress  increments  from  the  principal 

strain  increments  has  been  done  according  to  the  method  derived  by  Lubahn 
(41 

and  Felgar.v  '  This  appendix  outlines  the  method  and  gives  the  assumptions 
used  in  the  computation. 

The  following  assumptions  are  necessary  in  order  to  make  the 
computation: 


1) 

2) 

3) 


Acjg  *  0;  a  condition  of  plane  stress  pertains 
aa.  »  Ae^  +  Ae^  (i  =  1  or  2  for  the  principal  strains) 
The  deformation  theory  of  plasticity  pertains,  so  that, 
with  Aag  =  0 

Ae-j-Ae2  Ae2_4e2  Ae^-Ae^ 

Acr-j-At^  A0£  Aa-j 


(Al) 


(A2) 


4)  During  plastic  deformation,  plastic  strain  does  not 
cause  a  volume  change,  so  that 

Ae-j  +  AE2  +  Aej  3  0  (A3) 

5)  The  elastic  stress-strain  relations  are: 


Aa-j 

4e!  =  T  •  E  4o2 


Aao 

Aeg  =  - £-  Aa-j  (A4) 

Ae^  3  ^  (Aa-j  +  A02) 

These  equations  may  be  combined  and  reduced  to  a 


quadratic  equation  in  Ao^1 


(A5) 


a  Ac2  +  b  A0£  +  c  *  0 


where  a  ■  l-2u  b  a  -E ( 2ax -j  +  AAg ) 
so  that 


c  *  -[Aa.|  a  +  Aa-j  b] 


Ac.- 


.b-/b2-4ac 


‘2  2a 

The  computational  procedure  for  calculating  the  principal  stress 
increments  from  the  principal  strain  increments  is  as  follows: 


(A6) 


1 


2. 

3. 

4. 

5. 

6. 

7. 

8. 


If  both  ax i  and  AA2  are  less  than  the  elastic  limit  strain 
in  a  tensile  sample, Ac^  and  Aag  are  computed  from  Equation  (A4) 
For  all  other  cases,  steps  2-9  are  followed. 

Assume  a  value  of  Aa-j. 

Compute  Ao2  from  (A6). 

Compute  =  (Aa-j  -  Ac-jAf^  +  At^  ) 

Compute  Ae^  from  (A4). 

Compute  Ae..  from  (Al). 

2_ 

ST 


1/2 


(A7) 


1/2 


2  2  2  •  - 

Compute  (Ae-j  +  Ae-|Ae2  +  Ae2  )  .  (A8) 

Determine  Aa®^  from  the  appropriate  stress-strain  relation 
for  the  material*.  For  fatigue  cracks,  an  equation  of  the 
following  form  has  been  used: 

Ae  m 

Aaeff  “  K1  (A9) 


★ 

The  cyclic  stress-strain  relation  is  normally  obtained  from  a  specimen 
stressed  along  one  axis,  while  the  fatigue  crack  creates  a  multiaxial  stress 
state.  It  has  been  assumed  that  the  cyclic  stress-strain  curve  may  be  re¬ 
lated  to  the  multiaxial  case  through  use  of  the  effective  stress  and  strain. 
This  assumption  has  been  examined  by  Lambda  and  Si  debottom (5, 6)  who  found 
no  reason  to  Invalidate  it. 


9.  Compare  Aa*ff  with  If  Aa^f  -  Ao®ff  is  greater 

than  a  specified  tolerance,  assume  another  value  of  a<^  and 
return  to  step  2.  A  new  value  of  Aa-|  is  determined  from 

6 

-  Au  « 

AO'j  3  Aa-j  ( — “)  (A10) 

Aceff 
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8. 


SYMBOLS 


A  s  increment  of  a  quantity 

e.  *  elastic  strain  (i  =  x,  y,  z  or  1,  2,  3) 

e,j  *  plastic  strain 

3  total  strain 

eeff  *  effective  strain,  defined  by  Eq.  (A8) 

a.j  s  stress 

3  effective  stress,  defined  by  Eq.  (A7) 

E  3  Young's  modulus 

Poisson's  ratio 


u 


(a)  Minimum  load  level,  K  *  2.7  MN/m3;/2 


(b)  Maximum  load  level,  K  *  10.7  MN/m3^2 


Figure  1.  Photographs  of  the  crack  tip  for  which  the  analyses  in  this 
paper  are  made.  7075-T6  aluminum  alloy  grown  in  vacuum  at 
AK  =  8  MN/m3/2,  R  =  0.25. 
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I  . 


Figure  2.  Displacements  measured  from  the  photographs  shown  in  Figure  1. 

Notice  that  the  length  of  the  displacements  is  magnified  rela¬ 
tive  to  the  spatial  distances  shown.  The  dot  marks  the  unde¬ 
formed  location  and  the  end  of  the  line  displacement  due  to 
the  loading  increment. 
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DISPLACEMENT 


(b)  Fine  grid 


Figure  2  (continued).  Displacements  measured  from  the  photographs  shown 

in  Figure  1.  Notice  that  the  length  of  the  dis¬ 
placements  is  magnified  relative  to  the  spatial 
distances  shown.  The  dot  marks  the  undeformed 
location  and  the  end  of  the  line  displacement  due 
to  the  loading  increment. 
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Figure  4  (continued).  Maximum  values  of  the  stresses  computed  from  the  strains  of  Figure  3  in  MPa. 

Distances  along  the  x  and  y  axes  are  in  micrometers.  The  crack  is  shown 
schematically  on  the  plane  of  zero  stress. 
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